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ABSTRACT 

We present a high-resolution set of adiabatic binary galaxy cluster merger simulations using FLASH. 
These are the highest-resolution simulations to date of such mergers using an AMR grid-based code 
with Eulerian hydrodynamics. In this first paper in a series we investigate the effects of merging on 
the entropy of the hot intracluster gas, specifically with regard to the ability of merging to heat and 
disrupt cluster "cool-cores." We find, in line with recent works, that the effect of fluid instabilities 
that are well-resolved in grid-based codes is to significantly mix the gases of the two clusters and to 
significantly increase the entropy of the gas of the final merger remnant. This result is characteristic 
of mergers over a range of initial mass ratio and impact parameter. In line with this, we find that 
the kinetic energy associated with random motions is higher in our merger remnants which have high 
entropy floors, indicating the motions have efficiently mixed the gas and heated the cluster core with 
gas of initially high entropy. We examine the implications of this result for the maintenance of high 
entropy floors in the centers of galaxy clusters and the derivation of the properties of dark matter 
from the thermal properties of the X-ray emitting gas. 

Subject headings: galaxies: clusters: general — X-rays: galaxies: clusters — methods: N-body simu- 
lations, hydrodynamic simulations 



1. INTRODUCTION 

The growth of structure in the universe is thought to 
proceed in a bottom-up fashion, wit h smaller structure s 
merging into larger structures (e.g.. I White et al.lll993f ). 
The largest structures that have formed to date via this 
process are clusters of galaxies. Clusters of galaxies are 
important astrophysical objects, both in their own right 
and for their use as cosmological probes. Galaxy clus- 
ters contain a good representation of the different forms 
and kinds of matter within the universe. The bulk of the 
mass in clusters of galaxies is comprised of dark matter 
(Zwicky 1937; Bahcah & Sarazin 1977), which is believed 
to be largely collisionless. The bulk of the baryonic ma- 
terial in clusters of galaxies is comprised of a hot dif- 
fuse plasma c alled the intracluster medium (ICM) (e.g., 
iSarazinl l 1988V The remaining (and smallest) component 
of mass is that of the galaxies themselves. Galaxy clus- 
ters provide an opportunity to witness the interplay be- 
tween these different kinds and forms of matter in close 
quarters. Additionally, clusters have become very useful 
as probes of cosmology. The number density of clusters 
of a given mass at a given redshift is a sensitive function 
of the parameters of the standard big bang cosmological 
model. Accurate measurements of the density and tem- 
perature structure of galaxy clusters enable determina- 
tion of cluster masses (under the assumptions of spheri- 
cal symmetry and hydrostatic equilibrium) which in turn 
helps to c onstrain the values of cosmological parameters 
(jKitavam a & Suto 1996; Voit 2005). 

Extensive study of both of these aspects of astrophysics 
with galaxy clusters requires state-of-the-art observa- 
tions and simulations. The latest generation of X-ray 
telescopes {Chandra^ XMM-Newton) has provided high- 
resolution observations of the ICM that have made it pos- 
sible to probe sensitively the density, temperature, and 



metallicity structure of the gas. These observations have 
made possible the discovery of a myriad of features in the 
ICM such as shocks, cold fronts, and bubbles, to name a 
few. In addition, such precise measurements have made 
hydrostatic estimates of cluster masses a useful tool for 
constraini ng cosmological parameters on a par with other 
methods (jVikhlinin et al.l l2009b). In tandem with these 
advances in observation, astrophysical simulation codes 
have taken advantage of large-scale computing resources 
and innovative algorithms to model clusters of galaxies 
with unprecedented resolution and with the ability to in- 
clude detailed models of the relevant physical processes 
occurring in clusters. 

One of the most fundamental of these physical pro- 
cesses that shapes the life of clusters is merging. The 
merging process is responsible for forming the clusters 
and greatly influences the state of the matter within a 
cluster throughout its lifetime. Mergers drive shocks 
into the intracluster medium, heating and compress- 
ing the gas and driving turbulent motions. As a re- 
sult of merging the kinetic properties of the dark mat- 
ter are also significantly altered. Mergers also drive 
galaxy activity, producing new bursts of star forma- 
tion as the gas is compressed, or in s ome cases pos- 
sibly inhibiting it (jFuiita et al.l [1999"; "Be kki fc CoucE 
l2003f ). Many observed clusters of galaxies show signs 
of current or recent merging. Two prominent examples 
of rec ent note are IE 065 7-56 (the "Bull et Cluster") 
(Mark evitch et al.l 12 002: C lowe et al.l I2006D and Abeh 
520 (Mar kevitch et al.n2005l : iMahdavi et"aLll2007[ ). Ex- 
plaining the myriad of observable consequences of merg- 
ers using simulations is an active area of research. 

From the perspective of simulations, there are two ways 
to study the merging process. The first is by studying 
mergers that occur in the context of a cosmological simu- 
lation, "out in the wild." The advantage of this approach 
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is that it captures as accurately as possible the merging 
process as it occurs in the real universe, with many merg- 
ers of many different mass ratios and impact parameters 
going on at the same time, in the context of the over- 
all cosmological expansion. The second computational 
approach to studying mergers is by focusing on particu- 
lar merger scenarios in isolation via simulations of ideal- 
ized cluster mergers. These simulations involve the setup 
of two or more idealized, spherically symmetric clusters 
in a localized setting and with an initial orbital trajec- 
tory. Though this method is useful for studying individ- 
ual a ct ual merg:ers as seen in obs er vations (e.g;. iT akizawa 
^0061; 'SDrinsel & Farrar' "2007'; 'Akahor i Yoshikawa 
2008; MastroDietro & Burkert 2008; Rand all et alll2008l : 
ZuHone et al.. |2009a,^), an additional application has 
been to explore a space over parameters such as the mass 
ratio and initial impact parameter of the progenitor clus- 
ters or by varying the physics of the cluster ICM or dark 
matter. 

Several previ ous investigations of merger parame- 
ter s pa ces (e.g. iRoettiger et al.l 119971: [Ricker fc SarazinI 
2001"; "Gomez et al. 120021: i Ritchie &: Thomas 2002": 
McCaxthv et al. 200^; iPoole et al. 2006, 2007, 2008; 
Takizawa et al.l 120 lOt ) have emphasized the effects of 
merging on cluster morphology, the thermodynamics 
of the gas, and observables. We build on this work 
with a new set of simulations. The set of simula- 
tions presented here is the highest-resolution parame- 
ter study of idealized binary cluster mergers performed 
with an AMR-based PPM code to date, with a resolu- 
tion 3-4 times higher t han the most recent investigation 
(jRicker fc SarazinI I200L) based on such a formalism. 

In this first of two papers, we use these simula- 
tions to examine the effect of merging on the en- 
tropy of the cluster gas. Observational surveys of 
clusters have shown that although in nearly all clus- 
ters the entropy (defined here as 5* = ksTne ^^^) 
follows a power-law profile with r adius, the central 
core entropies of clusters v ary (Llovd-Davies et al.'"2000'; 
Pon man etlil l 120 031: rV^iFet al 2005: Faltenbacher et al 
2007; Morandi fc Ettoril 120071 : iCavagnolo et al.l 120091 : 
Pratt et al.. >2009.) . These observed entropy profiles are 
well-fit by an entropy profile of the form 

Where is the value for the co re entropy and 

a 21 1.0 — 1.3. In particular, Cavagnolo et al.l 

(|2009f ) and iPratt et al.l p009) showed that clusters ex- 
hibit a bimodal distribution of central core entropies. 
ICavagnolo et al.l (j2009,) found a low-entropy peak at 
^0 ~ 15 keV cm^, a high-entropy peak at 5*0 ~ 
150 keV cm^, and very few clusters with core entropies 
5'o - 30 - 50 keV cm^ Similarly, IPratt et al.l (|2009i) 
found peaks at 5*0 3 and 5*0 ~ 75 keV cm^ . These core 
entropies are higher than what would be expected from 
a simple ^^cooling flow" scenario (e.g. Fabian fc NulserJ 
[19771 : lFabianl[T99l . where the core gas is allowed to cool 
and compress unabated until it turns into stars or molec- 
ular clouds. Hi^h-resolution XMM-Newton spectroscopy 
(jPeterson et al." 2001", "2003) and Chandra spectral imag- 
ing (David et al...2001.) have shown that there is indeed 



little gas below T 1 keV in the cores of clusters with 
some of the highest cooling rates. 

Since cooling via X-ray radiation is directly observed, 
but the expected central entropies from unabated cool- 
ing is not, clusters with high entropy floors require 
a heating mechanism. Proposed candidates for heat- 
i ng of cluster cores inc lude magnetic fleld reconnection 
(jSoker fc Sarazinlll990f ). thermal conduction due to elec- 
tron collisions (e.g., Narayan fc Medvedev 2001; Fabian, 
Voigt, fc Morris 2002; Zakamska fc Narayan 2003) and 
turbulent conduction (e.g., Cho et al. 2003; Voigt fc 
Fabian 2004), and heating by cosmic rays (e.g., Inoue 
fc Sasaki 2001, Colafrancesco fc Marchegiani 200 8 ); a re - 
cent review can be found in iPeterson fc FabianI ()2006f ) . 
The currently favored mechanism is heating by the cen- 
tral AGN (e.g., Bohringer et al. 1993; Binney fc Tabor 
1995; McNamara et al. 2001, 2005; Fabian et al. 2006; 
Forman et al. 2007; for a recent review see McNamara 
fc Nulsen 2007). The AGN explosions blow the bub- 
bles in the ICM and inject energy into the ICM in the 
form of relativistic particles as well as mechanical energy. 
However, the precise mechanism by which this energy 
heats the central ICM is still unclear. A flne balance be- 
tween AGN explosions and cooling is required to avoid 
the complete blow-up of the cool cores, which gave rise 
to "feedback" models, where the cooling flow itself feeds 
the AGN. 

Thermal conduction is another particularly attractive 
idea, because it taps the vast reservoir of thermal energy 
in the gas just outside the cool core, while automatically 
ensuring that the core will not be overheated, since the 
heat influx decreases with diminishing temperature gra- 
dient. The classic plasma conductivity via Coulomb col- 
lisions was shown to be insufficient even at its full Spitzer 
value (e.g., Zakamska fc Narayan 2003). It has a strong 
temperature dependence and decreases right where it is 
most needed, and tangled magnetic fields should fur- 
ther suppress it (as was indeed observed outside the cool 
cores, Markevitch et al. 2003a). Additionally, magneto- 
hydrodynamic instabilities such as the heat-fiux-driven 
buoyancy instability (HBI) may suppress heat conduc- 
tion to cluster cores from the cluster outskirts alto- 
gether by aligning magnetic field lines perpen dicular to 
the ra dial temperature gradients in clusters ([Q uataert] 
2008; Parrish fc QuataertI 120081 : iBogdanovic et a l. 2009; 
Parris h et al.l l2009l ) , though recent works have shown 
that modest amounts of turbulence in cluster cores may 
work to keep the fields randomized an d sustain some con- 
duction of heat to the c luster core (jRuszkowski fc (jhl 
I20T0I : iParrish et ani20Tol ). Though the effects of AGN 
and conduction may account for the lack of systems with 
extremely low central cooling times and entropies, it is 
also unclear whether or not either effect can explain the 
highest entropy floors seen in some clusters, which can 
reach central entropy values of ~ 100 or higher. 

Previous investigations of cluster merge rs (e.g., 
Ritchie fc T homas 120021 : iMcCarthv et al.l 120071 : 
Poole et al.l |200 8) have indicated that so-called cluster 
"cool cores" (characterized by temperature inversions, 
high gas densities, and low central entropies and cooling 
times) are largely resilient to the heating that occurs as 
a result of merging. The central entropies and cooling 
times increase for a short time after the core passage, but 
the heating is not enough to offset cooling, which quickly 
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reest ab lishes a cool core within a few Gyr (iPoole et al.l 
I2QQ8D . iMcCarthv et all l\2007\ ), in a suite of adiabatic 
idealized merger simulations, determined that the ICM 
is heated in two major stages: the first is when the 
shocks are generated at the initial core passage and the 
second is at a later stage when material is shock-heated 
as it recollapses onto the merger remnant. Though the 
overall entropy of the merger remnant was increased 
with respect to the progenitor systems consistent with 
the self-similar prediction, the power-law behavior of 
the entropy profile was maintained all the way down 
into the central core (see Figure 5 of that paper). It is 
notable that these SPH merger simulations also resulted 
in little mixing of the cluster gases, since such mixing 
could provide a source of heating to cluster cores if the 
gases are of very different entropies. 

It is known that the various methods for solving the 
equations of hydrodynamics result in different degrees of 
mixing. Whil e previous investigations of mixing due to 
merging re.g.. lRitchie fc Thomasll2QQ2l : lPoole et al.ll2QQ8f ) 
have employed the Lagrangian "smoothed-particle hy- 
drodynamics" (SPH) formulation for solving the Euler 
equations, our simulations have been performed using an 
Eulerian PPM formulation. The two different approaches 
to hydrodynamics differ in their ability to appropriately 
mo del turbulence and the onset of fluid instabilities (see, 
e.g. iDolag et al.ll20Q5l : lAgertz et al.ll2QQ7l : IWadslev et al.l 
The amount of mixing of the gas is also dependent 
on the precise nature of the physics of the ICM, as the 
presence of magnetic fields, turbulence, and dissipative 
processes will affect the mixing of the gas. The extent 
to which the ICM mixes as a result of mergers and its 
effect on the thermodynamic properties of the cluster gas 
is therefore dependent the physical processes operating 
in the ICM as well as the chosen algorithm for solving 
the hydrodynamic equations. Therefore, it is important 
as a baseline to characterize the mixing process and the 
effects of such mixing on the thermodynamics of the clus- 
ter gas in the simplest model for ICM, that of an inviscid, 
unmagnetized gas, in the context of an Eulerian, PPM- 
based hydrodynamics code. Doing so is the goal of this 
paper. 

This paper is organized as follows. In Section [2] we 
describe the simulation method, setup, and rationale for 
the initial conditions. In Section [3] we present the re- 
sults of the simulations, focusing on the degree of mixing 
of the intracluster gas components and the alteration of 
the thermodynamic state of the cluster g as as a result. In 
Section m we discuss the implications of these results for 
different approaches to hydrodynamics for galaxy cluster 
simulations, the survival of cool cores in cluster mergers, 
and the use of X-ray observables to determine properties 
of the dark matter. In Section [5] we summarize this work 
and suggest avenues for future investigations. Through- 
out this paper we assume a spatially flat ACDM cosmol- 
ogy with h 0.7, 0.3, and Qb Omh''^. 

2. SIMULATIONS 
2.1. Method 

We performed our simulations using FLASH, a parallel 
hydrodynamics/7V-body astrophysical simulation code 
developed at the Center for Astrophysical Thermonu- 
clear Flashes at the University of Chicago (^Fryxell et all 



[2000h . FLASH uses adaptive mesh refinement (AMR), 
a technique that places higher resolution elements of the 
grid only where they are needed. In our case we are 
interested in capturing sharp ICM features like shocks 
and "cold fronts" accurately, as well as resolving the in- 
ner cores of the cluster dark matter halos. As such it 
is particularly important to be able to resolve the grid 
adequately in these regions. AMR allows us to do so 
without needing to have the whole grid at the same res- 
olution. FLASH solves the Euler equations of hydrody- 
namics usin g the Piecewise-Parabolic Method (PPM) of 
iColella fclVoodwardI (p,984), which is ideally suited for 
capturing shocks. FLASH also includes an A/'-body mod- 
ule which uses the particle-mesh method to solve for the 
forces on gravitating particles. The gravitational poten- 
tial is computed using a multigrid solver included with 
FLASH (|Rickerll2008[ ). 

2.2. Initial Conditions 

For these idealized cluster merger simulations, we 
choose initial conditions based on information gleaned 
from cosmo logical simulations and cluster observations. 
Since we are interested in the effects of merging on re- 
laxed systems, our initial clusters are configured to be 
systems consistent with observed relaxed clusters and 
cluster scaling relations. For the temperature normal- 
ization, we choose clusters that lie along the M500 — Tx 
relation of IVikhlinin et all ()2009af ). which is 



M500 = MoE{z)-\Tx/5 keV)" 



(2) 



with Mo = 3.02 x IO^^H-^Mq and a = 1.53. The depen- 
dence on the redshift z is encapsulated in the function 
E{z) = H{z)/Ho. For the total gas mass within r 5oo , we 
use the relation (also from IVikhlinin et al.ll2009ai) 

fg{h/0.72y-^ = 0.125 + 0.037 logioMis (3) 
where M15 is the cluster total mass, M500, in units of 

For our grid of simulations we have chosen a set of 
masses representative of galaxy clusters and groups. The 
chosen masses are given in terms of M2005 correspond- 
ing to the mass within which the average density is 
200pcrit(^)- The mass M200 is related to the radius r2oo 
via 

47r, 



M200 = --[200pcrit(^)]r: 



200 



(4) 



M200 is calculated from M500 from the assumed value of 
C200 and the assumed form of the density profile. 

For the total mass distribution (gas an d dark mat- 
ter) we choose an NEW (tNavarro. Frenk. fc White II1997D 
profile: 

^^^^^'^^ = / / (5) 

r/rs(l + r/rs) 

Where the scale density ps and scale radius are deter- 
mined by the following constraints: 



^200 
C200 
200 3 

2 ^200Pcrit(^j 



log(l + r/rs) 



r/rs 



1 + r/rs 



(6) 



(7) 
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Table 1 

Initial Cluster Parameters 



Cluster 


M200 


^200 


C200 






^0 


^1 






(Mq) 


(kpc) 






(keV) 


(keV cm2) 


(keV cm2) 


Np 


CI 


6 X 1014 


1552.25 


4.5 


0.1056 


4.97 


9.62 


192.40 


5,000,000 


C2 


2 X 1014 


1076.27 


4.7 


0.0879 


2.42 


5.08 


101.60 


1,684,119 


C3 


6 X IO13 


720.49 


5.1 


0.0686 


1.10 


2.73 


54.60 


513,137 



Table 2 

Initial Merger Parameters 



Simulation 


R 


b (kpc) 


V^200,l 


A 


SI 


1:1 


0.0 


0.0 


0.0 


S2 


1:1 


464.43 


0.3 


0.013 


S3 


1:1 


932.28 


0.6 


0.026 


S4 


1:3 


0.0 


0.0 


0.0 


S5 


1:3 


464.43 


0.3 


0.009 


S6 


1:3 


932.28 


0.6 


0.018 


S7 


1:10 


0.0 


0.0 


0.0 


S8 


1:10 


464.43 


0.3 


0.003 


S9 


1:10 


932.28 


0.6 


0.006 



The NFW functional form is carried out to r2oo. For 
radii r > r2oo the mass density fohows an exponential 
profile: 



Ptot(r) = 



Ps ^ V X ( ^ ~ ^^^^ ' 

C200(l + C200)^ V^200 y V ^decay , 



(8) 



where is set such that the density and its first derivative 
are continuous at r = r2oo and rdecay = 0.1r2oo. 

The total mass distribution uniquely determines the 
gravitational potential, and via the constraint of hydro- 
static equilibrium of the gas within this potential, the 
pressure profile is uniquely determined. To determine 
the physical quantities of density and temperature, more 
information is needed. Since we are here interested in 
the effect of mergers on the entropy profile of clusters, we 
make use of Equation [1] to set the gas properties in our 
initial systems. We start off with small core entropies 
(^o, see Table [1]) in order to make our models consis- 
tent with relaxed, "cool-core" galaxy clusters, and set 
a = 1.1. With this information in hand, we write the 
condition for hydrostatic equilibrium as follows: 



dP 

dr 

d{PgT) 



d 



jinip dr 



IJLvap dr 

3/2' 

t' ' 



TV 
S 



s) 1^ 



(9) 
(10) 

(11) 



Following other authors (jPoole et al.ll2QQ6[ ) , we solve this 
differential equation by imposing T(r2oo) = ^^200, where 

GM2ooprnp 



200 



2r2oo 



(12) 



is the so-called "virial temperature" of the cluster. With 
the temperature profile in hand we may use it and 
the entropy to determine the gas density profile, and 



so also determine the dark matter density profile via 

Pdm = Ptot — Pg- 

After the radial profiles are determined it remains to 
set up the distribution of positions and velocities for the 
dark ma tter particles. Here we fo llow the procedure out- 
lined in iKazantzidis et aP (|2QQ6l ) . For the particle po- 
sitions a random deviate u is uniformly sampled in the 
range [0,1] and the function u = MDM(^)/^DM(^max) is 
inverted to give the radius of the particle from the center 
of the halo. For the particle velocities, the procedure is 
less trivial. Many previous investigations have made use 
of the "local Maxwellian approximation." In this proce- 
dure at a given radius the particle velocity is drawn from 
a Maxwellian distribution with dispersion cr^(r), where 
the latter quan tity has been derived from solving the 
Jeans equation (jBinnev fc Tremainelll987l ). It has been 
shown that this approach is not sufficient to accurately 
represent the velocity distribution functions of dark mat- 
ter halos with a central cusp such as the NFW profile 
(Kazantzidis et al. 2004). To accurately realize particle 
velocities, we choose to directly calculate the distribution 
function via the Eddington formula (^Eddington . 1 9 1 6f ) : 



1 



\/87r2 



d-^c 



d^ 



_]_ / dp 

^^£ 



^=0 



(13) 

where ^ = — <l> is the relative potential and £ = ^ — ^v^ 
is the relative energy of the particle. We tabulate the 
function T in intervals of £ interpolate to solve for 
the distribution function at a given energy. Particle 
speeds are chosen from this distribution function using 
the acceptance-rejection method. Once particle radii and 
speeds are determined, positions and velocities are deter- 
mined by choosing random unit vectors in 3?^. The sim- 
ulation parameters for each simulated cluster are given 
in Table [H 

For all of the simulations, we set up the two clus- 
ters within a cubical computational domain of width 
L = 10/i~^ Mpc on a side. The distance between the clus- 
ter centers is g i ven by the sum of their respective r2oo. 
IVitvitska et al. (|2002l ) demonstrated from cosmological 
simulations that the average infall velocity for merging 
clusters is 'yin(^vir) = 1-1^7 where Vc = \/ GM^i^/r^i^ is 
the circular velocity at the virial radius rvir for the pri- 
mary cluster. For all of our simulations, this is chosen 
as the initial relative velocity {v^^ ~ 1200 km/s). The 
clusters are initialized in the x — y coordinate plane at z 
= 0. 

For this investigation we explore a parameter space 
in mass ratio of the clusters R and in initial impact 
parameter b. The set of simulations in this parameter 
space is detailed in Table O Taking the virial mass of 
our primary system as M200 = 6 x lO^^M©, we per- 



Figure 1. Schematic representation of initial merger geometry. 



form mergers involving virial mass ratios of R = 1:1, 
1:3, and 1:10. We al so va ry the impact parameter b 
for our runs. Vitvits ka et al.l (|2QQ 2) also demonstrated 
that for merger encounters the average tangential veloc- 
ity is vj_ ~ 0.4^0- Instead of splitting our initial veloc- 
ity into radial and tangential components, we choose to 
initialize our systems with impact parameters that re- 
sult in relative tangential velocities consistent with the 
results of Vitvitska et al. (2002). Since the velocity v± 
represents an average velocity, we select a range of tan- 
gential velocities consistent with this relation. Figure [1] 
shows the initial merger geometry and how the tangen- 
tial and radial components of the velocity are related to 
the initial choice of impact parameter. Finally, for each 
merger, we detail the angular momentum parameter A 
= J|^|1/VGM^/2 (jPeeblesI [19691) as a measure of the 
angular momentum of each simulation in Table [2l 

The highest AMR resolution of this set of simulations 
is Ax = 5h~^ kpc, which is ^3-4 times higher than the 
resolution in Ricker fc SarazinI (|200l[ ). The number of 
dark matter particles in each simulation is given in Table 

HI 

Since we will be concerned mainly with the effects of 
binary merging on cluster cores, it is important to ensure 
that the dramatic changes that will be seen to occur in 
the simulations are the result of merging and not evo- 
lution of the isolated primary cluster initial conditions. 



To this end we have run a single-cluster test, evolving 
the primary cluster in isolation for the same 10 Gyr as 
our merger simulations. Figure [2] shows the radial pro- 
files of gas density, gas temperature, dark matter den- 
sity, and gas entropy at the two epochs t = 0.0 Gyr and 
t = 10.0 Gyr. There is some moderate evolution in the 
profiles, particularly in the innermost few zones of the 
cluster core and in the cluster outskirts. The former is 
a consequence of force smoothing, (a known numerical 
effect due to the inability to resolve the gravitational 
force on scales smaller than the grid resolution). It is 
important to note that though there is some increase in 
the central entropy due to this effect, this increase is in- 
significant when compared to the increase that results 
from the mergers, as will be seen below. The evolution 
in the cluster outskirts is a consequence of the tapering 
off of the cluster profiles beyond r200 7 but has a very 
small effect on the profiles only in these regions and does 
not affect the cluster core. 

3. RESULTS 

3.1. Qualitative Description: Entropy Maps 

Over the course of a galaxy cluster merger, entropy 
is generated via processes that dissipate kinetic energy 
into heat energy, such as shocks and turbulence, and pro- 
cesses that bring into contact gases of different entropies 
and mix them. Here we show maps of the ICM entropy 



6 




throughout the merger to get a broader picture of the 
major entropy-generating events during cluster mergers. 
Each set of maps are taken from shces through the z = 
coordinate plane, the plane of the cluster centers' mutual 
orbit. 

3.1.1. Equal-Mass Mergers 

For our equal-mass mergers, there is no distinction be- 
tween the primary and the secondary system, so either 
cluster could be chosen as the primary without any loss 
of specificity (in all discussions that follow of the equal- 
mass mergers, it is assumed that the cluster whose center 
of mass has x < is the primary). In the head-on case 
(simulation SI, Figure [3]), the cluster cores approach each 
other, driving equal and opposite shocks into the gas. 
As the dark matter cores pass each other, the shocked 
gas in between the clusters forms a flattened low-entropy 
"pancake" -like structure (t ~ 1.6 Gyr). As the shocks de- 
part, this gas expands and cools adiabatically for a time, 
but eventually falls back toward the center of the sys- 
tem. The gas at high radii perpendicular to the merger 
axis develops Ray leigh- Taylor instabilities and is quickly 
mixed in with the surrounding gas. Meanwhile, the cen- 
tral gas is dragged back and forth by the oscillating dark 
matter cores {t ~ 3-4 Gyr). This process drives more 
shocks into the gas, heating it further. Eventually the 



gas forms a core with a roughly constant entropy of 6* ~ 
300 keV cm^ with a radius of r ~ 300 kpc. 

The off-center equal-mass cases (simulations S2 and 
S3, Figures [4] and [5|) are very similar to each other. In 
these cases the gas cores do not collide at the first core 
passage but sideswipe, creating a bridge of stripped, low- 
entropy gas that stretches between the two dark matter 
cores (t ~ 2 Gyr). The gas remaining in the cores is ini- 
tially held back by ram pressure as it moves outward, but 
as it moves out to lower density regions the effect of this 
pressure is weakened and the gas "slingshots" through 
and out the other side of the potential well. This gas 
expands adiabatically, curving around the dark matter 
cores and cooling, forming two large plumes of gas. The 
heads of these plumes fall back along with the dark mat- 
ter cores to the center of the system by t ~ 3 Gyr and are 
heated. The plumes themselves are stable for a few Gyr 
but eventually (by t ^ 4 Gyr) they have become mixed 
in with the surrounding gas due to the Kelvin- Helmholtz 
instability. The central gas in these cases forms a roughly 
constant-entropy core of 5* ~ 500 keV cm^ with a radius 
of r ~ 400 kpc. 

3.1.2. Unequal- Mass Mergers 

The unequal mass mergers have a qualitatively dif- 
ferent evolution, as a result of the asymmetry in mass. 
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Figure 3. Slices of entropy through the z 
1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 8.0, and 10.0 Gyr. 



10.0 56.2 316.2 1778.0 10000.0 

log S (keV cm^) 

= coordinate plane for simulation SI [K = 1:1, h 
Each panel is 5 Mpc on a side. 



kpc). The epochs shown are t = 0.0, 



which results in the formation of similar but somewhat 
different structures. In the head-on, 1:3 mass-ratio case 
(simulation S4, Figure [6j), the denser, lower-entropy core 
of the secondary cluster punches through the less dense 
core of the primary (t ~ 1.5 Gyr), disrupting it com- 
pletely. As the secondary's gas moves forward, a conical 
stream of low-entropy gas that has been stripped from 
the core trails it, falling into the gravitational potential 
well of the primary cluster {t ~ 2.0 Gyr). As this gas 
falls in, the coherent structures are ripped apart by fluid 
instabilities and mixed in with the surrounding, higher- 
entropy ICM. At the same time, the secondary's core. 



which had been initially trailing the dark matter core 
due to the ram pressure of the primary's gas, now en- 
counters regions of lower density. This causes the gas 
to fall back into the dark matter core, which then over- 
shoots the potential minimum and begins to climb out of 
it on the other side, expanding and cooling adiabatically 
{t ^ 2.0-3.0 Gyr). At this point, the secondary's dark 
matter core has already begun to fall back towards the 
primary's. By t = 4.0 Gyr the gas remaining from the 
secondary's core has formed a low-entropy, collimated 
stream that flows into the coalescing dark matter halos. 
By the end of the simulation at t = 10.0 Gyr this gas has 
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Figure 4. Slices of entropy through the z ■ 
as in [3] Each panel is 5 Mpc on a side. 



10.0 56.2 316.2 1778.0 10000.0 

log S (keV cm^) 

: coordinate plane for simulation S2 (R = 1:1, b = 464 kpc). The epochs shown are the same 



thoroughly mixed in with the higher-entropy gas of the 
primary and has formed a roughly constant entropy core 
with 5'o - 450 keV cm^ 

In the off-center, unequal-mass cases (simulations S5 
and S6, Figures [71 and [8]), a very different situation 
emerges. The secondary cluster does not disrupt the core 
of the primary cluster initially, but the gravitational dis- 
turbance caused by its passage sets up sloshing of the 
central core gas in the dark matter-dominated potential 
{t ^ 2.0 Gyr). This is due to the fact that as both the 
dark matter and the gas are pulled toward the secondary 
cluster during the first core passage, the gas is held back 



by ram pressure and becomes displaced from the poten- 
tial minimum. As the gas sloshes back into the potential 
minimum and then out again, it expands and cools, but 
as it comes into contact with higher entropy gas at larger 
radii it mixes with this gas and is heated. Meanwhile, 
the gas from the secondary cluster has formed a plume 
and stripped-gas bridge connecting back to the primary 
in much the same way as in the equal-mass cases. The 
bridge accretes cold gas back onto the primary cluster, 
followed by the plume which falls in which is mixed in 
with the sloshing central gas {t ^ 4-6 Gyr). Through- 
out this time period the plume is relatively stable but is 
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Figure 5. Slices of entropy through the z ■ 
as in [3] Each panel is 5 Mpc on a side. 



10.0 56.2 316.2 1778.0 10000.0 

log S (keV cm^) 

: coordinate plane for simulation S3 (R = 1:1, b = 932 kpc). The epochs shown are the same 



eventually broken up and mixed in with the surrounding 
gas. 

In our most extreme mass ratio cases, the effect of the 
subcluster on the main cluster is smallest, but the results 
are qualitatively similar to the 1:3 mass ratio cases, with 
the associated core disruption and accreting stream in 
the head-on case (simulation S7, Figure [9]) and subclus- 
ter plumes and main core sloshing in the off-center cases 
(simulations S8 and S9, Figures [TOl and [TT]). However, 
in these cases the onset of the effects of instabilities on 
these features is much sooner after the merger than in 
the previous cases. The Kelvin- Helmholtz instability be- 



gins to break up the plumes and sloshing core cold fronts 
almost immediately following the core passage. This is 
because the gravitational support from the subcluster is 
considerably weaker, making it much easier to strip the 
gas and providing less gravitational support against the 
onset of instability. 

3.2. ICM Mixing 

In order to investigate the degree of mixing of the gas of 
our initial systems, we have made use of "mass scalar s" . 
These quantities are defined as Yi = pi/ p, where pi is the 
partial density, and are advected along with the fluid: 
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Figure 6. Slices of entropy through the z ■ 
as in [3] Each panel is 5 Mpc on a side. 



10.0 56.2 316.2 1778.0 10000.0 

log S (keV cm^) 

- coordinate plane for simulation S4 (R = 1:3, 6 = kpc). The epochs shown are the same 



dt 



V • (pFiv) = 



(14) 



For the primary and secondary systems we define a 
mass scalar for each that is initially set equal to unity for 
radii within r2oo- We then examine the degree of mixing 
of these two components at later times. To parameter- 
ize this mixing we follow iRitchie fc ThomasI (|2002f ) and 
define the degree of mixing as 



M = 1 



Pi - P2 



Pi + P2 



(15) 



Given the above definition for the mass scalar s, 
relation simply reduces to 



this 



M = 1 



Y1-Y2 



(16) 



For example, a value of M = 1.0 indicates equal 
amounts of the two cluster gas components in a cell 
(Y1/I2 = 1), a value of M = 0.5 indicates that the clus- 
ter gas is mixed in a ratio of Y1/I2 = 3 or Y1/Y2 = 1/3, 
and a value of M = 0.0 indicates the gas is not mixed at 
all, and that only the gas from one cluster is present. 
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Figure 7. Slices of entropy through the z ■ 
as in [3] Each panel is 5 Mpc on a side. 



10.0 56.2 316.2 1778.0 10000.0 

log S (keV cm^) 

: coordinate plane for simulation S5 (R = 1:3, b = 464 kpc). The epochs shown are the same 



Figure [12] shows the resulting mixing M with gas den- 
sity contours overlaid at t = 10 Gyr, in a slice though the 
z = coordinate plane. Additionally, Figure [13] shows 
the mixing M as a function of radius. For the equal 
mass mergers, the gas of the two clusters is well- mixed 
in all cases of varying impact parameter out to a radius 
of nearly '^0.5 Mpc. In the head-on case (simulation 
SI), the gas is also significantly mixed along the plane 
perpendicular to the merger trajectory where the large 
"pancakes" of gas formed. The degree of mixing is more 
pronounced in the off-center cases (simulations S2 and 
S3). In these cases the fact that the cores slide pass each 



other and are stripped before finally merging into a single 
remnant allows for more mixing between the two than in 
the case where the collision is head-on. For the 1:3 and 
1:10 mass ratio mergers, the gas is also mixed but it be- 
comes less efficient as the impact parameter increases. In 
the head-on case with a 1:3 mass ratio (simulation S4), 
the disruption of the primary's core by the secondary al- 
lows for significant mixing of the two cluster components. 
The result is very different for the off-center cases, how- 
ever. In simulation S5 {b = 464 kpc) in the center of the 
cluster the two components are only mixed in a 3-to-l 
ratio and in simulation S6 {b = 932 kpc) the central gas 
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Figure 8. Slices of entropy through the z ■ 
as in [3] Each panel is 5 Mpc on a side. 



10.0 56.2 316.2 1778.0 10000.0 

log S (keV cm^) 

: coordinate plane for simulation S6 (R = 1:3, b = 932 kpc). The epochs shown are the same 



is hardly mixed with gas from the other cluster at all. 
Most of the mixing in these cases occurs at radii 
500 kpc, as the stripping of the secondary's gas occurred 
as it made its first core passage, and as this gas fell to- 
ward the center of the primary it was broken up by insta- 
bilities and mixed with the surrounding gas. Similarly, 
in the head-on case for the 1:10 mass ratio (simulation 
S7) the central gas is mixed in a 3-to-l ratio, but in the 
off-center cases (simulations S8 and S9) the central gas 
is not mixed with gas from the other cluster at all. In 
these extreme cases, the low entropy core is not disrupted 
completely but only displaced from the center somewhat 



and as the system as a whole relaxes this material sinks 
back to the center of the newly-reconstituted potential 
well. 

Figure [M] shows snapshots of volume-rendered 
movies of gas mixing in our galaxy clus- 
ter merger simulations that can be found at 
|http://www.cfa . harvard.edu / ~ j zuhone / Research /Entropy] 
_and_Mixing.html. ^ ^ 

3.3. Radial Profiles 

It is instructive to examine the final radial profiles of 
various quantities compared to the initial radial profiles 
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Figure 9. Slices of entropy through the z ■ 
as inO Each panel is 5 Mpc on a side. 



10.0 56.2 316.2 1778.0 10000.0 

log S (keV cm^) 

: coordinate plane for simulation S7 (R = 1:10, 6 = kpc). The epochs shown are the same 



for the same quantities of the primary cluster. Whereas 
our initial gas profiles represent typical cool-core clusters 
with peaked gas density profiles, a positive radial tem- 
perature gradient at low radii, and power-law entropy 
profiles almost all the way to the cluster center, the final 
profiles appear very different. 

We have constructed these radial profiles by first tak- 
ing the cluster potential minimum as the cluster center. 
By visual inspection we find that at the end of each sim- 
ulation this position corresponds well with the location 
of the gas and dark matter density peaks. Secondly, we 
divide the cluster into spherical radial annuli of width 



Ar = 10.5h~^ kpc, wide enough so that there are enough 
points/particles in each annulus to adequately sample 
each relevant quantity, but narrow enough to sample the 
radial range adequately. We have tested the robustness 
of our profiles by varying the width of our radial an- 
nuli by a factor of 2-3x, for which we find no significant 
changes in our profiles. The velocity mean and velocity 
dispersion are estimated from the sample mean and vari- 
ance of the cells for the gas and the particles for the dark 
matter. 

Figure [15] shows the final electron number density pro- 
file for each simulation, with the initial density profile 
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Figure 10. Slices of entropy through the z ■ 
same as in [3] Each panel is 5 Mpc on a side. 



10.0 56.2 316.2 1778.0 10000.0 

log S (keV cm^) 

coordinate plane for simulation S8 (R = 1:10, b = 464 kpc). The epochs shown are the 



of the primary shown for comparison. In ah cases, the 
central density is reduced by a factor '^3-10, and the 
gas mass is redistributed such that the resulting den- 
sity profile is considerably flattened towards the center 
as compared to the initial density profile. Such gas den- 
sity profiles are better fit b y the traditional /3-model 
(jCavaliere fc Fusco-Femianol [l976) than those of "cool- 
core" clusters, in which the gas density continues to rise 
towards the center. 

Figure [16] shows the final temperature profile for each 
simulation, with the initial temperature profile of the 
primary shown for comparison. The initial profile of the 



primary is typical of a "cool-core" cluster, featuring a 
temperature inversion with a negative radial gradient in 
the outskirts to a positive radial gradient at radii r ^ 
200 kpc. Each merger converts the temperature pro- 
file from this form to one where the temperature rises 
almost continuously towards the center within approxi- 
mately the same radial range. The differences are dra- 
matic: in the equal-mass merger cases (simulations Sl- 
S3) the central temperature is increased by nearly AT ^ 
12 keV, with the 1:3 mass-ratio cases increasing it by 
AT ^ 7-9 keV and the 1:10 mass-ratio cases increasing 
it by AT ^ 3-7 keV. 
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Figure 11. Slices of entropy through the z ■ 
same as in [3] Each panel is 5 Mpc on a side. 



10.0 56.2 316.2 1778.0 10000.0 

log S (keV cm^) 

coordinate plane for simulation S9 (R = 1:10, b = 932 kpc). The epochs shown are the 



Figure [TTl shows the final gas mass fraction (/gas(< 
r) = Mgas(< r)/Mtot(< r)) profiles for each simulation, 
with the initial gas fraction profile of the primary shown 
for comparison. Each merger considerably lowers the gas 
mass fraction in the central regions (by nearly a factor of 
4 or more in each case), corresponding to the fiattening 
of the gas profile, distributing it to lar ger radii. 

Recent works (e.g. Hansen 2009; Host et al. 2009; 
^altenbacher et al. 2007) have suggested that the combi- 
nation of X-ray observations with some simple assump- 
tions about the form of the dark matter velocity disper- 
sion profile might allow constraints to be put on quan- 



tities such as the magnitude of the dark matter velocity 
dispersion as well as the degree of its anisotropy. A sim- 
ple way of achieving this is by defining a "dark matter 
temperature" , 



(17) 



where (Jdm is the 3-D velocity dispersion of the dark 
matter particles. This quantity was first introduced by 
llkebe et all (|2004[ ) to derive the dark matter velocity dis- 
persion and compare it to the gas temperature. Figure 
[T8l shows the final dark matter temperature profiles com- 
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Figure 12. The degree of gas mixing (as defined in equation I16|) in each simulation, at a time t = 10 Gyr. At the top, from left to right, 
are SI, S2, and S3, in the middle, from left to right, are S4, S5, and S6, and at the bottom, from left to right, are S7, S8, and S9. Each 
panel is 4 Mpc on a side. 



pared to the initial profile. Unlike the gas temperature 
profiles, the shape of the profile is relatively unaffected by 
the merger though the overall normalization is increased. 

This temperature is not well-defined, due to the fact 
that there is no thermodynamic equilibrium for a col- 
lisionless gas (if it were valid, the Maxwellian assump- 
tion mentioned in Section [2^ for generating initial con- 
ditions would be accurate), but it may still be phe- 
nomenologically useful. Assuming that the velocity of 
an individual particle is determined solely by the grav- 
itational potential, the equivalence principle might lead 
one to expect a similar temperature profile between the 
dark matter and gas components of a galaxy cluster, e.g. 
^dm(^) = f^Tg{r). In this way the observed gas tem- 
perature may be used to determine the kinetic proper- 



ties of the dark matter. In the case of relaxed galaxy 
clusters, it may be sensible to assume ~ 1. Figure 
[T9l shows the temperature ratio radial profiles for each 
simulation. For our initial relaxed primary cluster, the 
equality Tgas = Tdm is valid to within 10% for most of 
the radial range, excepting the cluster core within about 
50-100 kpc. For the final merged clusters, over most of 
the radial range the assumption is still accurate to within 
10%, but within the inner 100-200 kpc the profiles begin 
to diverge. In the innermost regions of the final merger 
remnants, Tdm is only ~20-30% of Tgas- 

Our idealized progenitor clusters began with zero gas 
velocities relative to the rest frame of the cluster, so it 
is instructive to examine the velocities which are still 
present in our resultant systems at the end of our sim- 
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Figure 13. Gas mixing profiles at t = 10 Gyr. Left: Gas mixing profiles for the 1:1 mass-ratio simulations. Center: Gas mixing profiles 
for the 1:3 mass-ratio simulations. Right: Gas mixing profiles for the 1:10 mass-ratio simulations. 





Figure 14. The degree of gas mixing in 3-dimensional volume renderings, taken a t particular snapshots. Full movies of gas mixing can 
be found at |http://www. cfa.harvard.edu/ ^jzuhone/Research/Entropy_and-Mixing.htmH Left: Simulation S3. Right: Simulation S4. 



ulations. Figure [20] shows the total gas velocity disper- 



sion profile, ag 



<ji), for each simula- 



tion. Each merger converts a fraction of its energy into 
random motions of the gas, with the final velocity dis- 
persion ag ranging from 100-500 km s~^. By far, the 
head-on mergers generate more turbulent motion than 
the off-center mergers, which may have significant veloc- 
ity dispersions at larger radii r ~ 300 kpc but fail to 
penetrate the center of the primary, and the resulting 
dispersions are correspondingly less. Figure [2T] shows 
the mean radial velocity profile, ^7^, for each simulation. 
By the end of each simulation, the mean radial veloci- 
ties are quite small, with speeds \vr \ ^50 km s~l, with 
the largest velocities at the highest annuli of our profiles 
(r^l Mpc). Figure [22] shows the circular velocity pro- 
file at each radius, Vc = \f{vQ^ ^ V(^\ Not surprisingly, 
head on-mergers generate the least amount of circular 
velocity (14^100 km s~^), whereas off-center mergers 
induce more rotation, up to ~ 400 — 500 km s~^ 
in the strongest cases (equal mass mergers). The high- 
est circular speeds at the end of each simulation are at 
radii r ^ 100 kpc. Finally, Figure [23] shows the ratio of 
the velocity dispersion and the circular velocity to the 
sound speed of the gas squared; this provides a measure 
of the pressure support that may be provided by these 



motions. These values range from ~ 10-20% over the radii 
r < 200 kpc, with the highest levels of support reaching 
up to ~30-40% in the cluster outskirts. These levels of 
non-thermal pressure support from random and circular 
motions are consistent with previous results from cos- 
molog;ical simul ations (e.g., Piffaretti & Valdarnini 20081; 
[Fang et al.''2009l: ILau et al.> 2009). 

In our nine simulations we saw from the entropy maps 
that the final central entropy and the size of the entropy 
core varies as both the impact parameter and cluster 
mass ratio is varied. Figure [24] shows the final entropy 
profile for each simulation, with the initial entropy profile 
of the primary cluster shown for comparison. The com- 
bined effect of the heating and expansion of the central 
gas of the cluster is to raise the entropy of this central 
gas to a high adiabat. All of the final profiles exhibit a 
large entropy core where the entropy fiattens out in the 
center, out to a radius of r ~ 100 — 500 kpc, and the value 
of the central entropy is significantly larger than the ini- 
tial value, with ~ 100-500 keV cm^. The behavior of 
the profile outside of this central fioor is close to power- 
law, essentially unchanged from the progenitor clusters, 
consistent with the observational and simulation works 
that show that the real difference between clusters in this 
regard is in the height of the entropy floor. 

Finally, the entropy of the gas in the ICM is closely 






related to the cooling time, the total energy of the ^as di- 
vided by the energy loss rate. It is given bv (iZhang fc Wul 
[2QQ3h : 

— Kt)(w)*(to^)""<'« 

where Ue is the electron number density, ksT is the 



gas temperature in keV, and g is the frequency-averaged 
Gaunt factor. With some manipulation, the cooling time 
can be expressed as a function of entropy and electron 
density: 

U = 28.69 Gyr ( ( ^% lO-^lrn-^ ) ' 

(19) 
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Figure 18. DM temperature profiles at t = 10 Gyr, compared to the initial profile. Left: Temperature profiles for the 1:1 mass- 
ratio simulations. Center: Temperature profiles for the 1:3 mass-ratio simulations. Right: Temperature profiles for the 1:10 mass-ratio 
simulations. 
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Figure 19. Profiles of the ratio of the dark matter "temperature" to the gas temperature at f = 10 Gyr, compared to the initial profile. 
Left: Temperature ratio profiles for the 1:1 mass-ratio simulations. Center: Temperature ratio profiles for the 1:3 mass-ratio simulations. 
Right: Temperature ratio profiles for the 1:10 mass-ratio simulations. 



In Figure [25] we show the coohng time profiles for the 
nine simulations at t = 10 Gyr, with the age of the uni- 
verse marked on the plot for comparison. All of the cen- 
tral cooling times have been significantly increased, by a 
factor of ~6-20 over the initial value of ~1 Gyr at the 
very center of the cluster. From Equation [191 we see that 
this is consistent with the large increases in entropy in 
the cluster cores. 

4. DISCUSSION 

4.1. To Mix or Not to Mix? 

In previous investigations of ICM mixing using con- 
trolled merger simulations, it has been reported that 
the ICM component in these simulations does not effi- 
ciently mix. Ritc hie fc ThomasI (|2002|), using SPH sim- 
ulations, defined the degree of mixing as we have above 
and reported that the ICM only mixed well in the cen- 
ters of merger remnants. iPoole et all (|2QQ8l ) investi- 
gated mixing from the perspective of metallicity gradi- 
ents, and found no significant fiattening of the metal- 
licity profile (which would imply mixing of gas of dif- 
ferent metalli cities) following a merger in any of their 
simulations. iRicker fc SarazinI (|2QQ1[ ) did not investi- 
gate mixing, and in that case the reduced resolution of 
their simulations would have resulted in the suppression 
of smaller- wavelength unstable perturbations, lessening 



the amount of mixing. Additionally, we are concerned 
with the effects of merger heating on progenitor clusters 
with cool cores, whereas the initial systems in that study 
were more akin to our resultant systems. On the other 
hand, in an idealized setup of an infa l ling subcluster into 
a large cluster of galaxies, iTakizawal (|2QQ 5) showed that 
as the subcluster falls into the core of the main clus- 
ter, its gas is thoroughly mixed due to Ray leigh- Taylor 
and Kelvin-Helmholtz instabilities. The resolution of the 
grid (Ax = 2 kpc) and the values of the chosen cluster 
masses in their simulations are similar to the values in 
our simulations. 

Our simulations demonstrate that significant mixing 
can occur as a result of major merging. Figures [3] 
through [TTI demonstrate that fiuid instabilities are ubiq- 
uitous throughout the duration of the merger, mixing 
the gas thoroughly. The presumed reason for the dis- 
crepancy is the choice of method to perform t he sim ula- 
tions. The s imulations of [Ritchie fc ThomasI (|2QQ2[ ) and 
[Poole et al.[ (|2QQ8. ) were performed using SPH methods. 
A recent comparison of galaxy cluster merger simula- 
tions from a SPH code (GADGET-2) and a PPM code 
(FLASH) revealed that in the SPH simulations mixing of 
the two cluster components is inhibited, whereas in the 
PPM si mulations the c luster components are thoroughly 
mixed ([Mitchell et al.[[2QQ9) . It has been demonstrated 
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Figure 20. Profiles of the gas velocity dispersion at f = 10 Gyr. Left: Velocity dispersion profiles for the 1:1 mass-ratio simulations. 
Center: Velocity dispersion profiles for the 1:3 mass-ratio simulations. Right: Velocity dispersion profiles for the 1:10 mass-ratio simulations. 





Figure 21. Profiles of the gas radial velocity at f = 10 Gyr. Left: Radial velocity profiles for the 1:1 mass-ratio simulations. Center: 
Radial velocity profiles for the 1:3 mass-ratio simulations. Right: Radial velocity profiles for the 1:10 mass-ratio simulations. 




Figure 22. Profiles of the gas circular velocity at f = 10 Gyr. Left: Circular velocity profiles for the 1:1 mass-ratio simulations. Center: 
Circular velocity profiles for the 1:3 mass-ratio simulations. Right: Circular velocity profiles for the 1:10 mass-ratio simulations. 



recently by lAgertz et al.l (|2QQ7f ) that spurious pressure 
forces can result at density gradients in SPH simulations, 
resulting from overestimations of the gas density of par- 
ticles approaching a high-density region. As a result the 
growth of Ray leigh- Taylor and Kelvin-Helm h oltz i nsta- 
bilities is inhib i ted. I n addition, iDolag et al.l (|2QQ5f ) and 
IWadslev et al.l (|2QQ8[ ) showed that the artificial viscosity 
used in most SPH implementations acts to damp tur- 
bulent motions. Therefore, in a galaxy cluster merger 



simulation, both of these effects can inhibit the mixing 
of gas from the two cluster components. In contrast, 
PPM simulations such as ours produce the right results 
for this physical setup, noting the fact that in real clus- 
ters other physical mechanisms (such as viscosity and 
magnetic fields) may come into play. The results of our 
investigation over a parameter space of mergers were an- 
ticipated^y the grid-based hydrodynamical simulations 
of lTakizawal (j2QQ5l) . In their simplified setup of two rigid 
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Figure 23. Profiles of the ratio of the gas velocity to the local sound speed at f = 10 Gyr, for both the velocity dispersion (solid lines) and 
the circular velocity (dot-dashed lines). Left: Profiles for the 1:1 mass-ratio simulations. Center: Profiles for the 1:3 mass-ratio simulations. 
Right: Profiles for the 1:10 mass-ratio simulations. 
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Figure 24. Gas entropy profiles at f = 10 Gyr, compared to the initial profile. Left: Entropy profiles for the 1:1 mass-ratio simulations. 
Center: Entropy profiles for the 1:3 mass-ratio simulations. Right: Entropy profiles for the 1:10 mass-ratio simulations. 




gravitational potentials, they showed that the gas from 
a subcluster falling into a large galaxy cluster will be 
ram-pressure stripped and mixed into the main cluster 
gas due to the action of Kelvin-Helmholtz and Rayleigh- 
Taylor instabilities. 

Supporting this, we find a correspondence between the 
amount of mixing of the two cluster gas components and 
the velocity dispersion of the gas in the final merger rem- 



nants (Figures [T3l and [2Q|). Where higher velocity dis- 
persions have been generated in the cluster cores, there 
is a corresponding higher degree of mixing. Figures [3] 
through [11] showed the ubiquitous presence of instabili- 
ties as each merger proceeds, and it is these instabilities 
that seeded the random velocities that efficiently mix the 
gas. In corresponding SPH simulations the instabilities 
would be largely absent. 
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Figure 26. Examples of core gas sloshing at late times in two of the offset mergers. The quantity plotted is gas entropy with the same 
color scale as Figures [3l through [TT1 Each panel is 3 Mpc on a side. Left: Simulation S5. Right: Simulation S9. 



A caveat to these results is that despite the high resolu- 
tion of these simulations, perturbations with wavelengths 
smaller than the smallest zone size {Ax ~ bh~^ kpc) are 
not resolved. Similarly, the turbulent cascade is not re- 
solved below this length scale. It is possible that resolv- 
ing these modes would enable more mixing of the gas 
than is present in these simulations. Simulations with 
high er resolution or subgrid r nodels for turbulen ce (see, 
e.g. iScannapieco fc Briiggenl I2QQ8I : iMaier et all [2009), 
would be necessary to determine whether or not this is 
the case. 

4.2. Destroying Cool Cores 

In all of the cases investigated here, the initially com- 
pact cool cores are completely destroyed as a result 
of the merger, and replaced with low-density, high- 
temperature, high-entropy cores. This is especially true 
in the cases where significant mixing of the cluster gases 
have occurred, such as each of the head-on mergers and 
all of the equal- mass cases. It is clear from Figures [3] 
through [TT] that this entropy increase primarily occurs 
through the mixing of high and low-entropy gases caused 
by the onset of fluid instabilities. 

There is a correlation in our simulations between the 
height of the central entropy floor and the degree of mix- 
ing. In the equal-mass mergers, the inner ~0.5 Mpc is 
very mixed, but there is a higher degree of mixing at 
larger radii in the off-center cases than in the head-on 
case. This is because in the off-center cases the cluster 
cores slip past each other and are stripped by instabili- 
ties, allowing for mixing over a large surface area, in con- 
trast to the head-on case where the cores collide inelasti- 
cally with each other and mixing only initially occurs at 
the boundary between the two clusters. Correspondingly, 
the entropy floor in simulations S2 and S3 is somewhat 
higher (nearly by a factor of 2) than in simulation SI, 
and the entropy floor has a larger extent. The differ- 
ences in mixing are most evident in the 1:3 and the 1:10 
mass-ratio mergers. In cases of zero impact parameter, 
the smaller, denser cluster (the secondary) penetrates the 



gas of the larger and allows for signiflcant mixing of the 
central ICM of the two clusters, nearly 100% for the R 
= 1:3, 6 = kpc case (simulation S4) and approximately 
60% for the = 1:10, 6 = kpc case (simulation S7). 
It is in these cases where the greatest increase in core 
entropy is seen. However, in the off-center cases, most 
of the gas from the secondary is stripped as it proceeds 
along its trajectory, and so, though there is mixing in 
the cluster outskirts, the core of the primary is consider- 
ably less mixed with material from the other cluster. We 
still see a signiflcant amount of mixing between the two 
components in simulation S5, so its entropy floor is com- 
parable to that of S4, but the entropy floor of S6 is lower 
(by about a factor of two) where the core of the primary 
has been relatively unmixed with that of the secondary. 
In simulations S8 and S9 the inner ~250 kpc has been 
relatively untouched, and they have roughly equal low 
entropy floors when compared to the merger product of 
simulation S7, which has been disrupted by the head-on 
collision. 

It should be noted that our diagnostic for mixing, M, 
only indicates where material from the different clusters 
has mixed together. Even in the simulations where the 
central cores of the primary cluster remain relatively un- 
touched by the gas from the secondary, mixing still plays 
an important role. The merger event sets of a process of 
"sloshing" of the cente r cluster gas in the da rk matter 
potential (as in Ascasib ar fc Mar kevitchi [20061 ). bringing 
higher entropy gas (which is still from the main cluster) 
in from higher radii, which is then mixed with the lower 
entropy gas, raising the central regions to a higher adia- 
bat (see Figure [26] for examples of core gas sloshing from 
our simulations). An example from these simulations of 
this effect is seen in the 1:3 mass-ratio cases, where the 
entropy floors of simulations S4 and S5 have been raised 
to the same adiabat but the gas from the two clusters has 
mixed only ~2/3 as much in the latter case. Although 
the mixing between the two cluster components was not 
complete, the central gas in this case still mixed with 
higher entropy gas from its own component, raising the 
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entropy in the central regions. A detailed examination 
of this process in the cont ext of an AMR code such as 
FLASH was carried out bv lZuHone et alJ (|201Qf ). 

In general, it has been shown that the results of non- 
radiative cosmological simulations depend on the numer- 
ical scheme adopted for hydrodynamics. Mesh-based 
Eulerian codes (such as FLASH) systematically pro- 
duce higher entropy gas cores than do particle-based 
hydrodynamic approaches (see , e.g;., iFYenk et all 119991: 
lO^Shea et all [20051 : iVoit et al.l 120051 ). [Mitchell et af 
([20091 ) concluded that the extra gas mixing in the PPM 
simulations was responsible for this difference. In their 
simulation the time period of greatest difference in en- 
tropy generation between the PPM and SPH codes was 
the timescale of merger evolution during which large 
vortices and turbulent eddies were present in the PPM 
simulation, mixing the gas. Supporting this conclu- 
sion, works where the SPH algorithm has been modi- 
fied to more properly resolve turbule nt motions, insta- 
bilities, and niixing of the gas (such as [Dolag et al.[[2QQ5[ : 
[Kawata et al.[[2QQ9[ ) results in higher central entropy pro- 
files in the merger remnants. 

Mitchell et al. (2009) noted that the ability of mesh- 
based codes to produce high entropy floors may have 
implications for observations of cluster entropy profiles. 
They speculated that heating from recent merger events 
might be responsible for the observed clusters that do 
not possess cool, low entropy cores, but rather large 
high-entropy cores such as the final merger products seen 
in our simulations. Re cent works (e.g. Cavagnolo et al. 
[2009[ : [Pratt et al.[ [2009[ ) have clarified the existence of a 
bimodal distribution of cluster central entropies, with a 
large population of clusters possessing low-entropy cores 
(5*0 ^30 keV cm^) and a similarly large population pos- 
sessing high-entropy cores 5*0 ^ 50 keV cm^). The former 
set of clusters corresponds to cool-core systems much like 
our initial systems. The latter set of clusters may have 
had their core entropies raised by r ecent mer^in^, and 
correspond to our resultant systems. 'Pratt pool ) 
pointed out a correlation between central gas mass frac- 
tion and central entropy, noting that the high-entropy 
floor clusters correspondingly had low central gas mass 
fractions. This is exactly what happens in each of our 
merger cases. This work indicates that mergers of clus- 
ters with even small subclusters could significantly raise 
the central entropy floor and central cooling time. Previ- 
ous investigations using astrophysical codes where mix- 
ing was suppressed have indicated that merging was not 
able to sustain high central entropies and cooling times, 
but simulations with realistic mixing such as these sug- 
gest the opposite conclusion. 

To demonstrate this conclusively, however, this work 
must be extended to more physically realistic setups. 
Mixing in the real ICM may not be as efficient or even 
more efficient than in our simulations due to additional 
physical processes that may be operational in the ICM. A 
real physical viscosity would act to dissipate turbulence 
and suppress fluid instabilities. Magnetic fields would 
only allow movement of the plasma along the field lines 
and would similarly partially suppress instabilities and 
mixing. Conduction, if efficient, would smooth out tem- 
perature, density, and vel ocity gradients, also sup press- 
ing instabilities (see, e.g. [Vieser fc Hensled[2QQ7l ). On 



the other hand, a more accurate characterization of the 
turbulent cascade (by either performing simulations at 
higher resolution or by implementing a turbulent sub- 
grid model) would likely increase the amount of mixing. 
Finally, to make a direct determination of the efficiency 
of heating generated from mergers cooling must be self- 
consistently included in the simulations. This would be 
necessary for a more direct comparison t o the results 
of cos mological simulations. For example, [Burns et al.[ 
^2008"), using a cosmological grid-based simulation with 
cooling and supernova feedback, have shown that cool 
cores are resilient to cluster mergers. In order to com- 
pare more directly with this result it would be necessary 
to include cooling in our simulations. Despite the limi- 
tation s of the pre sent work, it reinforces the conclusion 
of Mit chell et al.[ p009) that merging remains a serious 
candidate for heating the ICM. 

4.3. Deriving Dark Matter Properties from Gas 

Depending on the degree to which dark matter can be 
understood to be "collisionless" (in the absence of other 
mechanisms such as self-interaction or annihilation), it is 
not expected to be subject to the fundamentally micro- 
scopic processes of heating and cooling that affects the 
plasma of the ICM. Consequently, though the form of the 
DM "temperature" radial profiles are remarkably stable 
and retain their essential form after a merger (Figure (TH]), 
the gas temperature profiles have been transformed from 
cool-core like profiles with positive temperature gradients 
into profiles with negative or flat temperature gradients 
after a merger (Figure [T6|) . It therefore can be expected 
that methods of determining dark matter properties by 
assuming equivalence between the dark matter "temper- 
ature" and the temperature of the gas will fail in the cen- 
ters (r^200 kpc) of galaxy clusters. Figure [T9l confirms 
this expectation. In one sense, this reflects difficulties al- 
ready known with using the thermodynamic state of the 
ICM to derive properties of clusters, in that the large 
differences in the central temperature and density struc- 
ture between clusters (due to nongravitational heating of 
the gas such as AGN) must be accounted for in order to 
determine their effect on deriving quantities under the as- 
sumption of hydrostatic equilibrium. For example, global 
X-ray spectral temperatures of clusters are typically de- 
rived by excising the central emission out to some radius 
(typically 0.15r5oo, 0.1r2oo)7 in order to achieve a higher 
degree of simil arity between temperature profiles (e.g., 
Vikhli nin et al . 2006, 2009a). However, we saw that as 
a result of merging the deviation of K: from unity at the 
end of the simulation can be larger than 20% even out 
to a radius of r ^ 250 kpc, depending on the strength 
of the merger. Attempts to derive the kinetic properties 
of dark matter from the properties of the X-ray-emitting 
gas must take this fact into account when applying this 
technique to clusters that show signs of recent merging. 

5. CONCLUSIONS 

We have performed a set of fiducial binary galaxy clus- 
ter merger simulations that explore a parameter space 
over initial mass ratio and impact parameter. In this 
first of a series of papers we focus on the mixing of the 
intracluster medium between the two clusters and the 
resulting effect on the thermodynamics of the gas of the 
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merger remnant. In contrast to previous merger param- 
eter space studies but in agreement with recent works 
focusing on mixing, we find that the gas of the merging 
clusters mixes very efficiently. The degree of mixing is 
dependent upon the mass ratio and impact parameter of 
the cluster merger, but all of the mergers investigated 
here demonstrate a significant amount of mixing in some 
regions of the final merger remnant. This is a result of 
ubiquitous presence of fluid instabilities that occur as the 
merger generates large shear flows and sharp density gra- 
dients, which are resolved in grid-based fluid codes such 
as the one used in this study but not in particle-based 
fluid codes used in many previous works. These insta- 
bilities result in the mixing of high and low-entropy gas 
and as a result the thermodynamic state of the primary 
cluster's core is dramatically affected. In each case the 
cluster has been changed from one with a high-density, 
low-entropy, low-temperature core with a short cooling 
time to a low-density, high-temperature, high-entropy 
core with a long cooling time. This is the case even 
for mergers with smaller subclusters. Further evidence 
that mixing is the culprit is indicated by the correlation 
between final merged systems with high central gas veloc- 
ity dispersions and high central entropies, as the greater 
random motions in these systems allow for more efficient 
mixing. This result may also pose a difficulty for at- 
tempts to derive the kinetic properties of the dark matter 
from the thermodynamic properties of the X-ray emitting 
gas, which rely on the hypothesis of the near-equality of 
the gas temperature and dark matter velocity dispersion 
profiles. Additionally, such efficiency in heating the clus- 
ter core may be a primary factor in the generation of 
a population of galaxy clusters with high entropy floors 
and cooling times. Further simulation works including 
more realistic physics for the ICM (on the one hand, dis- 
sipative effects, magnetic fields, which would suppress 
mixing, and on the other hand, a more accurate char- 
acterization of the turbulence via higher resolution or 
a turbulent subgrid model which would enhance it) are 
needed to conffim this hypothesis. In spite of this, the 
main result of this work is that merging is still a serious 
candidate for the generation of heat to balance against 
cooling in the cores of galaxy clusters. 
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